Microscopic theory of vibronic dynamics in linear polyenes 
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We propose a novel approach to calculate dynamical processes at ultrafast time scale in molecules 
in which vibrational and electronic motions are strongly mixed. The relevant electronic orbitals 
and their interactions are described by a Hubbard model, while electron-phonon interaction terms 
account for the bond length dependence of the hopping and the change in ionic radii with valence 
charge. The latter term plays a crucial role in the non-adiabatic internal conversion process of 
the molecule. The time resolved photoelectron spectra are in good qualitative agreement with 
experiments. 

PACS Numbers: 31.25.Qm, 31.50.Gh, 33.60.-q, 31.70.Hq 



Photoexcitation in polyatomic molecules leads to the 
rapid mixing of vibrational and electronic motions, in- 
ducing charge redistribution and energy flow in the 
molecule. This non-adiabatic coupling is essential in 
photochemical processes in photobiological processes, 
such as those involved in vision Q , and in molecular elec- 
tronics [||. 

The progress in the technology of laser pulses has trig- 
gered the development of a new generation of spectro- 
scopies devoted to the investigation of ultrafast phenom- 
ena. Beautiful experiments have been recently reported 
on the study of ultrafast non-adiabatic processes in iso- 
lated molecules Time-resolved photoelectron spec- 
troscopy (TRPS) has been used to follow the vibronic 
dynamics of all- trans 2,4,6,8 decatetraene ||. In this 
molecule, the electronic ground state is the singlet So 
(l 1 A ff ), while the first optically allowed excited state, 
S2 (l 1 B u ), has an energy higher than the first excited 
"dark" singlet Si (2 1 A ff ). A laser pulse (pump) pre- 
pares the molecule in a vibrationally hot wave packet 
involving the state S2. The packet evolves and its time 
evolution is probed by photoexciting an electron and an- 
alyzing the ensuing spectra at subsequent times. The 
population of the "dark" band, associated with the Si 
electronic state increases with time, as a consequence of 
a non-adiabatic internal conversion between vibrational 
and electronic excitations. Such a process lasts a few hun- 
dred femptoseconds, and is a manifestation of the failure 
of the adiabatic description. 

The aim of any theory in this field is to relate exper- 
imental spectroscopical data with the basic microscopic 
interactions of the system. This task is particularly dif- 
ficult in the present case. A theory describing the exper- 
iments should start from a model which includes both 
electron-electron correlations as well as electron-phonon 
couplings, and should treat the latter beyond the adia- 
batic approximation. So far, to the best of our knowl- 
edge, time dependent photoemission spectra for such 
non-adiabatic internal conversion processes have never 



been calculated on the basis of a miscroscopic model. 
Previous theoretical work on the internal conversion pro- 
cess is based on semiempirical models for the relevant 
energy surfaces, supplemented by phenomenological cou- 
plings between them [Q. 

The goal of this Letter is to describe the whole process 
behind a TRPS experiment Q (laser pump followed by 
a dynamical internal conversion, probed by the ensuing 
time-dependent photoemission spectra), using a minimal 
microscopic model in which both electron-electron and 
electron-phonon interactions are exactly taken into ac- 
count. Our results show a good qualitative agreement 
with the experimental findings j^], and provide a link 
between these and the basic underlying microscopic in- 
teractions. 

The model we study is based on a formal separation of 
the hybridized s and p orbitals lying in the plane of the 
molecule (see Fig. 1(a)), and the p-7r orbitals perpendic- 
ular to it. The former are those with the larger contribu- 
tion to the chemical bond and vibration dynamics. The 
latter are those involved in the low-energy excitations of 
interest to us, and are described by a Hubbard Hamilto- 
nian. Its key ingredient is a strong local Coulomb inter- 
action U . Electron-phonon coupling terms are included, 
as derived from an expansion of the electron-ion interac- 
tion up to first order in the ionic displacements relative 
to the equilibrium positions |^|. The model reads: 
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n ia = fj CT /i,CT is the corresponding number operator, and 
qi = 1 — nj-f — Ujj is the net charge at site i. u\ are 
the ionic displacements with respect to the equilibrium 
positions, and, for simplicity, the vibration dynamics has 
been assumed to be that of a one dimensional chain of 
ions (the CH units) of mass m with a nearest-neighbor 
harmonic constant K. Two electron-phonon interaction 
terms have been included in Eq. ([!]), with coupling con- 
stants labeled by g and g' . The g-term is quite standard 
and accounts for the fact that the magnitude of electronic 
nearest-neighbor hopping becomes weaker or stronger as 
the bond stretches or compresses, respectively. The g'- 
term, with g' > 0, accounts for the contraction (expan- 
sion) of the bond upon removal (addition) of an electron. 
Such a term is a consequence of the dependence upon the 
valence charge of the spread of the wave functions (par- 
ticularly the a orbitals), a spread which in turn modifies 
the bond lengths. 

For sake of simplicity, and in order to be able to nu- 
merically solve exactly the Hamiltonian, we consider a 
molecule with L = 4 C atoms. The vibrational motion of 
the molecule for g — g' = can be described in terms of 
three normal modes, e , et, and e c , oscillating with fre- 
quencies uj a = J~(2 + V2)K/m, cub = \J {2 — V2)K/m, 

and oj c = ^/2K/m, as sketched in Fig. 1(b). The mode 
e;, will not be included, since it describes a uniform dilata- 
tion or contraction of the molecule, which should have a 
small coupling with the electrons in a larger molecule. 

In order to show that the Hamiltonian (0) represents 
the minimal model containing the relevant ingredients to 
describe the physics of decatetraene, we begin by dis- 
cussing the role played by each interaction term within 
the adiabatic picture. 

Role of U . The role of electronic correlations in the ex- 
citation spectrum of polyenes has been already discussed 
||,[l0|. Here if U were neglected, the optically allowed 
S2 state would have a lower energy than the dark Si 
state, while (U/t )>1.7 leads to the observed ordering of 
both singlet excitations. Given the reasonable estimate 
of to = 2 eV for the electronic hopping JsJ , we find that 
the value of U that best fits the data on decatetraene is 
U = 4.2 eV. 

Role of g. The electron-phonon interaction g is a stan- 
dard ingredient of the Pariser-Parr-Pople-Pierls (PPPP) 
model |9pC[|, crucial in explaining dimerization effects. 
The coupling between the electrons and the e a mode (the 
dimerization mode) leads to different equilibrium lengths 
for C-C and C=C bonds, alternating short and long 
bonds. From experimental data fill] we took uj a = 0.2 
eV (which implies K = 36 eV/A^and u c = 0.15 eV). 
A value for g — 3 eV/A is then calculated by fitting the 
experimentally observed difference in bond-lengths |pT| . 
Setting for the time being g' = 0, one can study the 
Born-Oppenheimer (BO) surfaces corresponding to the 
three lowest energy singlets So, Si, S2, as functions of 



the normal mode coordinates e a and e c (see Fig. 1(c)). 
The position of the minimum of the ground state BO sur- 
face, occurring at a finite e a > 0, reflects the dimeriza- 
tion. It is also evident from Fig. 1(c) that the electron- 
phonon coupling introduces anharmonicities in the BO 
energy surfaces leading to the occurrence of level cross- 
ings in the excited states, known as conical intersections 
. The S2 singlet is odd under space inversion while the 
Si is even. The conical intersection between both states 
is a central feature of the physics we want to describe. 
Note that the e a mode is even under inversion, while e c 
is odd. Then, on general symmetry grounds, the cou- 
pling between the electrons and the e c mode would be 
expected to produce a nonvanishing matrix element cou- 
pling the Si and the S2 states, leading eventually to an 
avoided crossing between the corresponding bands, with 
a gap between the two proportional to the effective inter- 
band tunneling amplitude ||. As a consequence of that 
mixing, the adiabatic description in terms of uncoupled 
BO surfaces would loose meaning, and the quantum na- 
ture of the phonons should be included explicitly. How- 
ever, in absence of the (/-term, the Hamiltonian possesses 
a subtle electronic symmetry (particle-hole P p -h) which 
leads to the vanishing of the interband tunneling matrix 
element, even for e c 7^ 0, as explained below. P p -h is 
defined as the invariance of Ti. (when g' — 0, and up to 
an inessential chemical potential shift) under the trans- 
formation /J — > (— iyfi t - a . Taking into account that 
the electronic states in polyenes correspond to half-filled 
configurations (the number of 7r-electrons N — L) and 
that this subspace is left invariant by P p -h, the eigen- 
states of Tt within this sector can be classified as even 
or odd according to P p -h B- Detailed analysis reveals 
that Si is even under P p ~h, while S2 is odd. As a conse- 
quence, the interaction term g alone cannot produce the 
interband coupling leading to the non-adiabatic effects 
observed in these molecules. 

Role of g' . For this reason, the additional interaction 
g' , which breaks particle-hole symmetry, is essential for 
the internal conversion in our model. Fig. 1(d) shows 
the BO energy surfaces when g' is included. Notice that 
the ground state energy surface is practically unaffected 
by this interaction. This is due to the fact that the on- 
site Coulomb repulsion inhibits charge fluctuations in the 
ground state, thus making the g' term, involving the net 
charges qi, ineffective. On the contrary, charge fluctua- 
tions are important in the excited states, and couple to 
the phonon modes through the g' term: As a consequence 
of the quantum tunneling for e c 7^ 0, due to g 1 , the two 
relevant potential-energy surfaces show now an avoided 
crossing (Fig. 1(d)), instead of intersecting (Fig. 1(c)). 
A rough estimate for the strength of g' can be obtained 
from its effect on a C2 dimer. By fitting experimental 
data for the bond lengths of C 2 (1.2425 A) and of 
(1.2682 A) with our model, we get g' » 3.95 eV/A [£JJ. 
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An alternative rough estimate based on the ionic radii 
of C+ 4 and C~ 4 ions (0.15 and 2.60 A respectively @), 
leads to g' w 10 eV/A. We found that s'«5 eV/A leads 
to reasonable values for the effective interband coupling 
and for the position of the conical intersection. 

In order to describe the observed internal conversion 
effects, the adiabatic picture must be abandoned. The 
normal coordinates are quantized as: e a = y // Ti/2mui a (a+ 
<r), e c — yj / h/2muj c (c + cr). The resulting Hamiltonian 
Ti can be numerically diagonalized by introducing a cut- 
off in the number of phonons n a and n c . We found that 
keeping states with up to n a — n c — 10 was enough to 
obtain accurate results. 

The laser pump is simulated by acting at t — on the 
ground state |^>o) with the following operator: 

6 = Y J G m \4> m ){(t> m \E, (2) 

m 

where \4> m ) denote the exact eigenstates of 7i, and the 
effect of the laser electric field is represented by the 
dipole operator E = ^ i Riqi 1 where Ri is the position 
of the C atom i for Ui=0. The pump pulse is assumed 
to have a Gaussian shape in time. The Fourier trans- 
form of its envelope determines the excitation ampli- 
tudes G m ~ exp (-(e m - e) 2 /2cr^), where e m = E m - E 
are the excitation energies relative to the ground state 
Eq, while e and a p are the mean excitation energy and 
the energy spread of the pump pulse. The Schrodinger 
time evolution of the prepared wave packet is then given 
by: \ip(t)) — exp(iHt/h)0\(j)o). A qualitative picture of 
the resulting dynamics can be obtained from P (t) = 
(t/j(t)\P \i/j(t)), where P Q is defined as the projector of 
the eigenstates of H on the subspace of electronic states 
with odd parity under space inversion. P (t) is shown 
in Fig. 2. The first two singlet excitation energies (with 
predominantly Si and S2 character) are ~ 3.74 eV and 
~ 4.31 eV respectively. The spectral density profile of 
the initial wave packet is shown in the inset of Fig. 2 
for a pulse with a p = 2 eV and e = 4.3 eV. The results 
for P (t) in Fig. 2 show that, on top of a fast oscillatory 
component, there is a slower internal conversion compo- 
nent which lasts approximately 243/s. This process is 
dominated by the two excited states labelled by (1) and 
(2) in the inset of Fig. 2. The excitation energies of these 
states are e\ = 4.304 eV and e 2 — 4.315 eV, respec- 
tively. State (1) has predominantly Si (even) character 
(((f>i\P \(f>i) ~ 0.4), while state (2) has predominantly S2 
(odd) character ((<f> 2 \Po\4>2) ~ 0.74). The evolution of 
a wave packet composed by just these two states is indi- 
cated with a dashed line in Fig. 2. This behavior suggests 
that the relevant time scale in the evolution of \ip(t)) is 
set by ~ h/(e 2 - d)- 

Within the sudden approximation, the photoelectron 
spectrum at time t is essentially a measure of the follow- 
ing spectral function |l5| 



pM = E E K#J/.>W*)>l a % - e'J, (3) 

where E' m and \4>' m ) denote the eigenenergies and eigen- 
states of Tt upon removal of a valence electron. Results 
for the time evolution of pt(e) are shown in Fig. 3, where 
an artificial broadening of the delta-functions has been 
introduced. The most salient feature of Fig. 3 is the 
transfer of spectral weight, as a function of time, from 
a group of states close in energy to the ground state of 
the system with N = L — 1 electrons (all of them with a 
predominantly Do nature, i.e. a doublet obtained by pho- 
toexciting S 2 ||), to a group of states at a higher energy. 
The latter states are identified as vibrationally hot states 
with predominantly Di nature (obtained by photoexcit- 
ing Si). The arrow in the lowest panel of Fig. 3 indicates 
the position of the lowest-energy state related to such 
a band. The behavior of the photoelectron spectra as 
a function of time is in qualitative agreement with the 
experimental results of TRPS. The larger energy gap be- 
tween the relevant photoelectron spectral features which 
we find in our simulation (~ 2 eV in our case, to be com- 
pared with the experimental result of ~ 1.2 eV ||) should 
be likely ascribed to the smaller size of the molecule we 
are considering in our simulation. 

In summary, we have shown that the model Hamilto- 
nian in Eq. ([[]) contains all the basic ingredients which are 
necessary to explain the essential features of the the dy- 
namics of linear polyenes, notably electron-electron and 
electron-phonon interactions treated in a non-adiabatic 
(fully quantum) framework. With such a model, we 
have described in detail the various features of a time- 
resolved photoemission spectroscopy experiment, obtain- 
ing a good picture of the underlying physics. We re- 
mark that, although the ground state properties as well 
as the structure of the electronic excited states are prop- 
erly described by the standard PPPP model ||[l0]], the 
proper treatment of the observed non-adiabatic internal 
conversion, related to tunneling between coupled Born- 
Oppenheimer states, requires the inclusion of usually ne- 
glected electron-phonon interactions, such as the g' term 
we have considered, which break a residual particle-hole 
symmetry of the Hamiltonian. On the technical side, we 
remark that exact diagonalization techniques are useful 
theoretical tools for the study of these effects, as shown 
by the reasonable agreement of our simulations with the 
behavior reported in Ref. j^]. 
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FIG. 1. (a) Scheme of the linear polyene. Carbon and 
Hydrogen atoms are denoted by C and H respectively, (b) 
Scheme of the normal modes of the molecule, (c) Cuts of the 
potential energy surfaces for g' = 0, with to — 2 eV, U = 4.2 
cV, g = 3 eV/A, K = 36 eV/A 2 . Dashed and solid lines 
correspond to e c = 0,0.1, respectively, (c) Same as (b) for 
g' = 5 cV/A. 
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FIG. 2. Evolution of the projector P a (t) (Solid line). The 
model parameters are those of Fig. 1(c). Dashed line: evo- 
lution of P (t) for a wave packet consisting of the two states 
(1) and (2) indicated in the inset. Inset: Spectrum of the 
prepared wave packet for a Gaussian laser pump pulse with 
e = 4.3 eV and o v = 2 eV. 
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